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Abstract 

Processes of electron-positron annihilation into charged pions and kaons are 
considered. Radiative corrections are taken into account exactly in the first order 
and within the leading logarithmic approximation in higher orders. A combined 
approach for accounting exact calculations and electron structure functions is used. 
An accuracy of the calculation can be estimated about 0.2%. 
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1 Introduction 

The processes in collisions of e + e~ beams at moderate high energies with detection of the 
final particles moving in centre-of-mass system (c.m.s.) at large angles are the subject of 
close attention at meson factories, such as VEPP-2M (Novosibirsk) M, DA$NE (Fras- 
cati) 0, BEPC/BES (Bejing) [||. The processes of pure quantum electrodynamics (QED) 
nature provide an important background for studies of subtle mesons properties. Besides, 
they may be used for a calibration and monitoring. Because of large cross-sections of 
the lowest order processes, radiative corrections (RC) to them are to be included in the 
consideration. 
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In our previous paper [[| we had considered the QED processes e + e~ — > e + e~( r y), 
/i + //~(7), 77(7). We developed there an approach for precise accounting of radiative 
corrections to differential distributions. The part of RC, describing the emission of hard 
additional photons was presented in the form convenient for imposing experimental con- 
ditions of the final particle detection. The contribution of higher orders of perturbation 
theory was considered in the leading logarithmical approximation. That was done by 
means of the structure function approach, writing a cross-section in the Drell-Yan form. 

In this paper we apply the same approach to processes e + e~ — > 7r + 7r~(7), KlKs(j), 
K + K^( r y), considering the pseudoscalar mesons as point-like objects. The effects of 
strong interaction of hadrons in the final state are parametrized by introducing form 
factors which are to be measured in an experiment. We assume as usually that vacuum 
polarization corrections (by hadrons and leptons) are also included in the form factors. 

The virtual and soft real photon emission corrections are calculated in the 0(a) order 
exactly. This permits us to keep explicitly the leading (containing large logarithm L = 
ln(s/mg), s = 4e 2 is the squared total energy in c.m.s.) and next-to-leading terms. 
The latters are regarded further as /C-factor terms in the Drell-Yan representation for a 
cross-section. 

Considering hard photon emission we extract contributions, containing large logarithm 
L. We distinguish the collinear kinematics of hard photon emission and the semi-collinear 
ones, which do not give rise to L. An auxiliary parameter, a small polar angle #0 1 with 
respect to the direction of the initial beams, is introduced for this purpose. The terms 
containing photon softness parameter A = Ae/s (Ae is the maximal energy of a soft 
(in c.m.s.) undetected photon) and the ones containing #0 from the contribution of the 
collinear region are regarded as compensating terms. When they being summed with the 
contribution of the semi-collinear region (with the relevant restrictions imposed) provide 
their finiteness in the limit A, 80 — > 0. 

This paper is organized as follows. In the second Section we consider the process of 
charged pion production. The explicit formulae (we keep pion mass exactly) for the lowest 
order virtual and soft real photon emission are presented. The charge-even and charge- 
odd contributions are given separately. The latter quantity permits us to obtain the charge 
asymmetry which can be measured. We put also the explicit formula for the differential 
cross-section of hard photon emission. In Sect. 3 a similar consideration is given for 
the case of neutral and charged kaon production near the threshold. In Conclusions we 
discuss the formulae obtained and their precision. The results are illustrated numerically 
in Figures. In the Appendix we discuss the structure of form factors and give the explicit 
expressions for the vacuum polarization operator. 

2 Pion Pair Production 

In the Born approximation the differential cross-section of the process 

e+(p+) + e-(p_) -> 7r + (g + ) + 7r(g_) (1) 

has the form 

^(*) = ^f *n 2 9\FM\', P=y/l- <M (2) 
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(p+ + P-) 2 = 4e 2 
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p q 



The pion form factor F n (s) takes into account vertex virtual corrections due to strong 
interactions and vacuum polarization by leptons and hadrons (including vector-meson 
resonances) 0] (see Appendix). We would like to underline that in our approach QED 
corrections are not included into F^s). One has to take the form factor from an experi- 
ment after an extraction of QED radiative corrections. 

Calculating QED radiative corrections we will consider pion as a point-like particle. 
We distinguish form-factor-type one-loop Feynman diagrams and the box-type ones. 
QED form factors are taken as F^ D (s) = 1 + F^ + 0(a 2 ). Using the known first order 
contributions to the electron and pion QED form factors 
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and the contribution due to soft photon emission |J [7], we obtain for the charge-even 
part of the differential cross-section the following formula: 
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Box-type diagrams and the interference of soft photon emission from electrons and 
pions give rise for the charge-odd contribution 
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The charge asymmetry in a quasi-elastic case has the form 

_ dq(c) - da(-c) _ da s J d v ^ 
dcr(c) + dcr(— c) dcr 

In the ultra-relativistic case {(3 — > 1) we obtain 
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This expression coincides with the result of Brown and Mikaelian || . 

The matrix element of the process accompanied by hard photon emission 



e (p_) +e + (p+) — 
can be presented in the following form: 
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(10) 



p_ + p+ = q- + q + + k, 



(p_ + p+) 2 , si = (g_ + g+) 2 , x± = P± k - 



Tensor T^ a describes the transition of a heavy photon into the system of two real pions 
and a real photon: 



7*(g) -> 7r + (g + ) + 7T (g_)+7(fe), g 2 = g 2 = m 2 



9 



£; 2 = 0. (11) 



Regarding CPT and gauge invariance the tensor can be written in the general form 
T; CT = a x m + a 2 m + a 3 L$ + A^CL, ^ = 0, IfTL = 0. (12) 



The last term (~ k a ) is irrelevant here. Tensors L^> read 

L ( pa = 9 P ° - k p q a , Lf} = qk Q p Q a - kQ (q a Q p + Q a k p ) + (kQ) 2 g pa , 

L fa = kQ ( y q 2 g pa -q p q a ) + Q a ( y qkq p -q 2 k p ), Q = ^(gi - q 2 ). (13) 

For the case of charge pion production, which is considered below, we used the approxi- 
mation of point-like pions, where 
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In reality some vector-meson resonance intermediate states give rise of contributions to 
the tensor. 

The differential cross-section of the process with hard photon emission reads 
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k± = k-p± , Q = ~(q + - q_) , 

P+P- 2 
t = -2p_q_, tx = -2p + q + , u = -2p_q + , u\ = -2p+q~, 

dr = d ^qlkP ^ -Q--Q+- k), 

where we use the notation: 

(abed) = -Sp abed = (ab) (cd) + (ad) (be) — (ac) (bd). (16) 
After algebraic transformations one can get a more compact expression: 
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In the ultra-relativistic limit (s ^> m 2 ) one has to drop in the above formula the terms 

^SS) ^SlSn ^SSi- 

For numerical estimations it is useful to separate the most singular part of the diffe- 
rential cross-section and to integrate it analytically. We mean the contribution due to 
hard collinear photon emission by electrons. We suggest the following procedure. Let us 
define narrow cones surrounding the momenta of the initial particles. Their opening angle 
is defined by an auxiliary parameter 6q. The vertex is taken in the interaction point. A 
photon emitted by the electron or the positron inside the cones (kp_ < 8q or kp + < 9q) 
will be called as a collinear one. On the parameter 8 one has to impose the restrictions 

1>0 O >— ■ (18) 

£ 

Integrating inside the cones we drop all terms proportional to 9 2 0. After simple calcu- 
lations one comes to a formula, where a factorization of a shifted Born cross-section can 
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be found. So, the process of collinear photon emission is factorized with respect to the 
hard process of the annihilation into pions. This is the manifestation of the known factor- 
ization theorem ||10|| . Indeed, we obtain the factorization of large logarithm L = ln(s/m 2 ) 
with an accompaniment of several non-leading terms. The dependence on the auxiliary 
parameter 6$ should cancel in the sum with the contribution of the integration outside 
the cones. We will keep ln(8 2 /A) terms and use them below as a compensator. 
The shifted cross-section da(zi,z 2 ) reads 
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where Z\ and z 2 are energy fractions of the almost real electron and positron after the 
emission of collinear photons. The energy fractions Yx 2 of the final pions can be found 
from the following kinematical relations: 
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The leading contributions to cross-section, containing large logarithm L, as may be 
recognized, combine into the kernel of Altarelli-Parisi-Lipatov evolution equation: 
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This formula is valid in the leading logarithmical approximation. We will modify it by 
including nonleading contributions and using the smoothed exponentiated representation 



for structure functions [11 
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In comparison with the corresponding formula in Ref. [TT] we shifted the terms, arising 
due to virtual e + e~ pair production corrections, form V 1 into T> e+e . 
The final expression for the corrected cross-section reads as follows: 
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Quantities a, 6 and k are defined in Eqs.(||,||). Schematically this expression can be written 



as 



da = der^ + [da™ + C^] + [da® + C®\ + da 
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where quantities C*- 2,3 - 1 denote compensators; and terms dex^ 2 ' 3,4 ) denote integrals of Ri,2,3, 
respectively. Starting with Eq.fll7D we can replace in the above expression quantities Ri 
in the following way: 
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Experimental conditions of the final particle detection are encoded by the 9-functions 
&cut(zi,z 2 ) (for the two-particle final state kinematics) and 9^ (for the three-particle 
one). They can be imposed explicitly by introducing the restriction of the following kind: 
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(25) 
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where y^e = e t h is the threshold of the detectors, angle ^ determines the dead cones 
around beam axes unattainable for detection. Angles Q± define the polar angles of the 
pions. More detailed cuts can be implemented in a Monte Carlo program W2 , using the 
formulae given above. 

There is a peculiar feature in the spectrum of hard photons. Namely in the end of the 
spectrum the differential cross-section is proportional to the factor 
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which defines a certain peak. It comes from the Feynman diagrams describing the emission 
by the initial particles ||13|| . 



3 Kaon Pair Production Near Threshold 



In the case of K^Kg meson pair production the differential cross-section in the Born 
approximation reads 

^M = <^L^e\F hS (s)\\ (27) 

Here, as well as in the case of pions production, we suggest that the form factor Fls 
includes also the vacuum polarization operator of the virtual photon. Quantity (3k = 



1 — 4m K fs is the i^-meson velocity in the centre-of-mass frame, and 9 is the angle 
between the directions of motion of the long living kaon and the initial electron. 
The corrected cross-sections has the form: 
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3, a ir 2 V 

1 + -6+ 

4 7T I 3 2 



24 V3 



-L - 2tt' 



37 
T 



- Ml 



+ 



+ 2 L * 



X 

~ 2, 

2 {- 

"2 1- 
.3 



1 



4(2 -x) In - + -(1 + 3(1 



In- 



x 



2m e 



x — 



(1 



(2 - 2x + x 2 



In 



sx 



it if 



1 — X 



+ {2-x) ln(l 



+ 



_ 5 N 
e(ar 



— 6 + x 

2 



b ( sx 
- m 



(29) 



mi 



1 
3 



2m f 



We omitted a small contribution (proportional to a{m$ — 2mK) / 'm^) from photon emission 
by the final particles. Note that at higher energies this effect is extremely interesting: it 
may shed light on the neutral kaons polarizability problem. 
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In the case of K + K mesons production the Coulomb final state interaction is to be 
taken into account: 

dQr 0( g ) a2 PK .2n\ p M |2 Z /onX 

z = *™ /s-4mf, / +8 -4m| r \- 1 

V V s \ si 



where w is the relative velocity of kaons |L4]. When s = m| we have u ~ 0.5. Because we 
consider the energy range close to <f> mass one may choose the maximal energy of the soft 
photon as 

y\ 1 

uj < AE = ms — 2mK *C %, A = « — . (31) 

25 

If required, more precise formulae for charge-even and charge-odd parts of cross-section 
may be obtained from the ones for charged pions production process with the replacement 



4 Conclusions 

Thus, we presented differential cross-sections to be integrated in concrete experimental 
conditions. The formulae are good as for semi-analytical integration, as well as for the 
creation of a Monte Carlo generator ||12|| . The idea of our approach was to separate 
the contributions due to 2 — > 2 like processes and 2-^3 like ones. The compensating 
terms allow us to eliminate the dependence on auxiliary parameters in both contributions 
separately. In the same approach we had considered in paper f| the processes of large- 
angle Bhabha scattering and electron-positron annihilation into muons and photons. 

Note that all presented formulae are valid only for large angle processes. Indeed, in the 
region of very small angles 9 « m e /e of final particles with respect to the beam directions 
there are contributions of double logarithmic approximation ||. These small angle regions 
give the main part of the total cross-section. We suppose that this kinematics is rejected 
by experimental cuts. 

Numerical computations were done for the case of pion production. In Fig.l we give 
the differential cross-section in the Born approximation, in the leading logarithmic ap- 
proximation (with the /C-factor) according to the first term of Eq. fl22|) , and as the total 
corrected cross-section. We used the following set of parameters: the energy threshold for 
pion registration Ai = y t h = 0.5, the detector angular acceptance 10° < 6± < 170°, the 
beam c.m.s. energy e = 0.51 GeV. The second Figure shows the compensation of auxil- 
iary parameters. Value of a® is shown by the solid line. It has to be summed with 
which is drown by circles. Value of 

i s shown by the dashed line. It has to be summed 
with drown by crosses. The parameters for Fig.2 are: Ai = 0.5, A = 0.01, 9q = 0.01, 
other parameters as in Fig.l. The corrected differential cross-section for charged pion 
production is shown in Fig.3. The beam energy and angular acceptance are the same as 
in Fig.l. For these numerical illustrations we did not take into account the pion form 
factor. 
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The precision of our results is defined by the omitted contributions. As concerns pure 
QED terms, the corresponding uncertainty is defined by unknown coefficients before the 
following terms: 



- L^ICT 4 , - ^10" 5 , e -L 2 ^10~ 7 , 32 

TV J \TV / TV S 

We estimate the unknown coefficient^ and derive the theoretical uncertainty of our calcula- 
tions for these three processes to be 0.2%. Some additional uncertainty is due to precision 
of the form factor determination in an experiment. There are also some model-dependent 
contributions due to different hadronic intermediate states, which are far beyond the scope 
of this publication. 
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Appendix 

Pion and kaon form factors are supposed to have the form i^(s) = (1— II(s)) .Fj(s), where 
the factors J i ' 7T ,ls,k{s) (unknown theoretically) take into account the strong interactions 
between hadrons in the final state. We put here the expressions for leptonic and hadronic 
contributions into vacuum polarization operator II(s): 
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(34) 



The first order leptonic contribution is well known [1]: 



ni(a) = \l - h - + f(x,) + f(x T ) 



- + 0(x M )6(l - x M ) + 0(x T )O(l - x r ) 



(35) 



: For small-angle Bhabha scattering we had performed the complete calculations of next-to-leading 
logarithmic terms []l6[ . The corresponding coefficients were found to be about unity. 
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where 
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^ + |(2 + x)v^^ln(i±^) for x < 1, 
| + |(2 + x)y/l - x atari (^7=1) for x > 1, 



(j){x) = \(2+x)^/l - X, X^ T = —. 

s 

In the second order it is enough to take only the logarithmic term from the electron 
contribution 

n 2 ( S ) = ~(L-i7r)+C(3)-A. (36) 

Here we present also a theoretical estimate for the contribution to the vacuum po- 
larization operator due to <fi meson in the energy region close to the resonance \fs ~ 

« 1020 GeV (see also || |17|). In this region one can write the cross-section of e + e~ 
annihilation into hadrons as follows: 

12* B„T% 

° k[a) = ( , - mir + m ' {) 

where B ee is the branching ratio of the decay — > e + e~, quantity is the total width 
of the meson. For the hadronic part of vacuum polarization connected with cf) meson we 
obtain: 

n fc (s) = n*(s) + n h ,(s), (38) 



where ELy(s) includes other hadronic contributions fla| . The contribution n^(s) is es- 
sential only in the region m§ — nT^ < */s < + nT^, n ~ 1. It has the following 
form: 

I± 

Ms) = gg= m \ B ee « 3.09 • 10- 4 . (39) 
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Figure 1: The differential cross-section da e e n /dc [nbarn] as a function of c. The 
Born cross-section is given by crosses. The cross-section in the leading logarithmic appro- 
ximation (with the /C-factor) is shown by circles. The solid line is the resulting corrected 
cross-section. 
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Figure 2: An illustration of the cancellation of the dependence on A and 8 . The solid line 
is a( 2 \ the dashed one is a^ 3 \ compensator C*- 2 -* is represented by circles and compensator 
(7(3) j s drown by crosses. 
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Figure 3: The value of the corrected cross-section da e e ->7r n /dc [nbarn] as a function 
of c and Ai. 
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